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ABSTRACT 



Areas of convective precipitation were delineated 
and averaged convective precipitation amounts determined 
within these areas. A statistical, step-wise, screening, 
linear regression procedure was used to correlate certain 
large scale parameters with -the convective precipitation 
amounts . 

Of the parameters considered in the two cases 
investigated, the best predictor of convective precipita- 
tion was either the terrain-induced vertical motion at 
the top of the boundary layer, or the product of the 
latter and the specific humidity at the same level* 



LIBRARY 

NAVAL P0SiC.R wJATS SCHOOL 
KONTERSY , CALIF. °C04Q 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 

ur\MT cdpv PA QA^AR-RIOI 



TABLE OP CONTENTS 



Section Page 

List of Tables 5 

List of Illustrations 7 

Acknowledgment 8 

1. Introduction 9 

2. Procedure 11 

3. Results of case studies 15 

4. Conclusions and remarks 29 

5. Bibliography 31 



3 



* 



LIST OP TABLES 



Table Page 

1. Matrix showing the simple correlation 21 

coefficients between the variables 

for Case A. 

2. Matrix showing the simple correlation 21 

coefficients between the variables 

for Case B. 

3. Results of step-wise regression analysis 22 

for Case A. 

4. Results of step-wise regression analysis 23 

for Case B. 

5 . Regression equations using only one 24 

variable . 

6. Matrix showing the simple correlation 24 

coefficients between the variables, 
computed for the beginning of the time 
period, for Case A. 

7. Matrix showing the simple correlation 25 

coefficients between the variables, 
computed for the beginning of the time 
period, for Case B. 

8. Regression equations using only one 25 

variable, which was computed for the 
beginning of the time period. 

9. Improved regression equations using 26 

only one variable. 



5 



LIST OF ILLUSTRATIONS 



Figure Page 

1. Surface synoptic situation for 1200 ig 

GMT 11 February, 1965 . 

2. Surface synoptic situation for 0000 19 

GMT 12 February, 1965 . 

3 . Surface synoptic situation for 1200 20 

GMT 16 March, 1965 . 

4. Surface synoptic situation for 0000 20 

GMT 17 March, 1965 . 

5 . Observed convective precipitation for 27 

Case A. 

6. Computed convective precipitation for 27 

Case A. 

7. Observed convective precipitation for 28 

Case B. 

8. Computed convective precipitation for 28 

Case Bo 



7 



ACKNOWLEDGMENT 



My gratitude is extended to Professor M. B 
for numerous helpful suggestions and criticisms 
the course of this research, and whose work was 
basis of this study. 



Danard 

during 

the 



8 



1 . INTRODUCTION 



The failure of numerical models to account for 
convective precipitation is well known (see, e.g., 

Danard (19 66 a, 1966 b ) ) . Such a shortcoming is expected 
since the equations used for numerical prediction are not 
applicable to sub-synoptic scale motion. Convective 
precipitation is highly variable in both time and space. 
This would prevent such precipitation from being accu- 
rately accounted for even with a much denser observational 
network than the one available in the United States today. 

This investigation is an attempt to delineate the 
areas of convective precipitation by the use of semi- 
empirical methods. Convective precipitation will be 
interpreted here as the average amount associated with 
sub -synoptic scale systems over a square area whose side 
is equal to the distance between adjacent grid points 
used in numerical prediction (381 kilometers at 60 degrees 
north). The amount defined is assumed to be determined 
by the large-scale properties of the atmosphere. 

Convective activity is confined to areas of unstable, 
or conditionally unstable, lapse rates. That is, the lapse 
rate should be steeper than the moist adiabatic lapse 
rate. Kuo (1965) describes a convective cloud model which 
will "bring surface air to all levels up to a great 
height so that inside the cloud the vertical distribution 
of temperature and mixing ratio are those of the moist 
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adiabat through the appropriate condensation level. 

This suggests that the rate of production of cloud and 
total water content would depend on the specific humidity, 
or mixing ratio, and the vertical velocity at some lower 
boundary. Ooyama (1963) assumed that the rate of 
precipitation was proportional to the flux of water 
vapor through the top of the frictional boundary layer. 

The total gain of water vapor in an air column, assuming 
no flux through the top of the column, may be represented 
by 

1- - (^1^/3 5 

(e„g., Krishnamurti (1968)), where I is the water gained, 
q is the specific humidity, and p is pressure. The sub- 
scripts b and t indicate the bottom and top of the column, 
respectively. Here the first term represents the flux of 
water vapor through the sides of a column of unit area, 
and the second term represents the flux through its base, 
which is taken here as the top of the planetary boundary 
layer. Since specific humidity decreases rapidly with 
height, the second term is frequently larger than the 
first if the level p^ is sufficiently high. 
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2 . PROCEDURE 



For the two case studies included here, Danard 
(1966b) has made numerical predictions of gridpoint 
precipitation within an area shown in Figures 5~8. 

Predicted amounts obtained from the first 12 hours of 
these integrations were assumed to represent precipitation 
effects of large-scale motions. Observed precipitation 
amounts were obtained by averaging all individual values 
reported by synoptic stations over a square centered 
on the grid point and having a side equal to the distance 
between adjacent grid points (approximately 381 kilometers) 
Convective precipitation amounts used are the difference 
between the averaged observed amounts and the amounts 
predicted by Danard, Only those points where the lapse 
rate was close to, or greater than, the pseudo-adiabatic 
value (see below) were used. The maximum 12 -hour 
precipitation amount predicted by the numerical scheme 
in these areas was only 0.4 millimeters in either case 
study. Thus, the areas of convective precipitation were 
clearly delineated from those of large-scale precipitation. 

A statistical, step-wise, linear regression procedure 
was used to correlate convective precipitation with the 
following parameters: (l) the difference between the pseudo 
adiabatic and the computed lapse rates (C"D> (2) the 
vertical velocity ( co ) at the top of the boundary layer, 

( 3 ) specific humidity (q) at the top of the boundary 
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layer, (4) the product of (2) and ( 3 ), and ( 5 ) dew point 
depression (T-T d ) in the lower troposphere. The program 
used is available in the Naval Postgraduate School 
computer library, where it is listed as BMD02R. A 
sequence of multiple linear regression equations is 
computed in a step-wise manner. At each step the 
variable is added which ma<kes the greatest reduction 
in the error sum of squares. 

For the two cases considered, average temperatures 
and dew-point temperatures were computed for the 850 , 

700, and 500 mb levels. The reported values at the 
beginning and end of the period were used to compute 
time averages. These averaged temperatures were used 
in the computations. 

All averaged parameters were computed from values 
reported by the stations. The resulting fields were then 
analyzed, and the values for the grid points used were 
taken from the analysis. 

As used here, the lapse rate is defined to be the 
difference between the temperatures at 850 and 500 mb. 

The pseudo-adiabatic lapse rate between 850 and 500 mb 
was taken from an Arowagram as a function of the tempera- 
ture at 700 mb. Convective precipitation was considered 
possible where the difference between the lapse rates 
»■ (r w ~D was less than 2.0°C/(350 mb). The areas delineated 
by this criterion included all regions where large amounts 
of convective precipitation was observed in both cases. 
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The dew point depression used was averaged between 
850 and 700 mb. Here the 500 mb level was omitted 
because, in both cases many stations reported motor- 
boating at that level. In the cases where motor- 
boating was reported at the lower levels,, the minimum 
dew point depression for which a radiosonde normally 
reports mot or -boating was used. 

Danard (1966a) expressed the vertical velocity at 
the top of the boundary layer as 

CO — - h 3 ( 2 ) 

where oj r and u)f are the vertical velocity induced by 
orographic and eddy-stress influences, respectively. 

In the area considered in this study (see Figure l) 
most of the terrain heights are less than 400 meters 
above sea level. Therefore, the surface geostrophic 
wind was taken to be that at 1000 mb. Thus, a>t is well 
approximated by 

ooo*^7Bf ) ( 3 ) 

where p^. is the standard atmosphere pressure at the 
elevation of the smoothed terrain provided by Berkofsky 
and Bertoni (1955 )> and ( ^g)]_Q00 t ^ ie S eostro Phi c 
wind at 1000 mb. 
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The frictionally induced vertical velocity was 



computed by 

v 3 ) • (4) 



Here V is a representative geostrophic wind in the 
planetary boundary layer and is the drag coefficient 
(Cressman (i960)). In the cases considered here ( V - )^qqq 
was used. However, over higher terrain, or in cases 
where ( Vg) 100 Q cannot be easily computed, ( V g )g^ 0 
might be more useful. 

Specific humidity at the top of the boundary layer 
may be expressed as 




o.6>ze e 

p - 0.3 78 e 3 




where e is the vapor pressure and p is the pressure in 
mib at the top of the boundary layer, assumed to be 85O 
mb here. Vapor pressures used were the saturation 
values (Smithsonian Institution ( 1951 )) corresponding 
to the dew point temperatures at 850 mb. 



>• 
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3. RESULTS OP CASE STUDIES 

Two synoptic situations were investigated; 1200 GMT 
11 February to 0000 GMT 12 February, 1965* and 1200 GMT 
16 March to 0000 GMT 17 March, 1965* These intervals will 
be referred to as Cases A and B, respectively. The 
surface synoptic situation for the beginning and end of 
the two 12 hour periods may be seen in Figures 1-4. 

Both cases were characterized by deepening cyclones 
moving from the southwest. These storms produced some 
of the worst weather of the 1965 winter in the central 
United States. 

The number of grid points used in Cases A and B 
(i.e., the sample sizes) was 16 and 17, respectively. 
Tables 1 and 2 give the simple correlation coefficients 
between all the variables used for Cases A and B, 
respectively. In both cases the correlation between a) 
and o)q was very high. In Case A, o)q had the highest 
simple correlation with convective precipitation (P). 

For Case B, the correlation between u) and P was slightly 
higher than that between u)q and P. The second highest 
simple correlation for Case A was between (T-T^) and P. 
However, there was little or no correlation between these 
variables in Case B. 

The results of the step-wise regression analysis may 
be seen in Tables 3 and 4 for Cases A and B, respectively. 
In the step-wise regression procedure, the first variable 



15 



entered is always the one with the highest simple 
correlation with P, unless the "deletion option" of 
the step-wise regression program is invoked. In Case 
A, after o)q entered the regression equation, u) was 
essentially redundant and did not enter the equation 
at all. For the same reason o)q did not enter the 
equation in Case B, since cJ entered first. In doth 
cases there was little improvement in the multiple 
correlation between the predicted value and P after 
the first variable entered the regression equation. 

An attempt was made to determine if a time delay 
between the flux of moisture through the boundary 
layer and resulting precipitation could be detected. 

To accomplish this, the regression procedure was 
repeated using variables computed for the beginning, 
rather than for the average, of the time period. The 
resulting simple correlation coefficients are given in 
Table 6 and 7* In Case A, there was noticeable improve- 
ment in the correlation between both oJq and P, and 
and P. However for Case B there was little change in 
these correlations. 

The step-wise regression program used permits the 
variable c*>q to enter the regression equation first for 
Case B if the variable u) is deleted. In this way it is 
, . possible to obtain the variable o)q as the primary 

predictor for both Cases A and B. In a similar manner, u) 
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has been selected as the primary predictor for Case A. 

The resulting equations, with only one variable entered, 
may be seen in Table 5- Due to the high simple corre- 
lation between a) and u)q there is little loss in explained 
variance for either case, when or uJq is employed as 
the primary predictor. 

In order to improve the prediction equations and 
to ensure that P^-0, predictions were made for both 
cases using the equation 



P = A 0 + A j (u)q) , f o r u < - Ao/Aj 

= O ) for ooq 'Z.-Aq/Ai 



The coefficients A q and A^ were determined by computing 
the quanity 

S 2 - I(P„ b -P) 2 , ( 7 ) 

where the summation extends over all the points in the 

sample. P ob is the observed convective precipitation at 

a gridpoint and P is the amount given by Equation (6). 

2 

S is simply the error sum of squares. Various values 

2 

of A q and A^ were tried until the minimum in S was 
found. This procedure was repeated using cO as the 
predictor. The improved prediction equations are shown 
in Table 9. Comparing the last columns of Tables 5 and 
9 , it is seen that the above procedure results in a slight 
reduction in error for Case B. However, for Case A the 
improvement is negligible. 
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Predictions were made for both cases using the 



equation 



P = «~ 0.45 - 1 .ZO(u)q) 5 for o)(\ < - 0.375 
= 0 j -for 0)3 Z — 0.375 




Equation (8) was obtained by averaging the coefficients 
given in the first two lines of Table 9 - The results 
may be compared with the observed convective precipitation 
in Figures 5-8. 
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Figure 1 
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Figure 3 




Figure 4 
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TABLE I 



Matrix showing the simple correlation 
coefficients between the variables for 
Case A, P is convective precipitation. 



p 


r w -r 


60 


q 


oOq 


T-T. 

a 


P 1.00 


0.28 


-0.43 


0.12 


-0.31 


-0.47 


r w - r 


1.00 


-0.34 


-0.13 


-0.35 


-0.44 


a) 




1.00 


-0.10 


0.95 


0.55 


q 






1.00 


-0.28 


0.01 


u)q 








1.00 


0.53 



TABLE II 

Same as Table I for Case B 





P 


r„-r 


CO 


q 


a)q 


T-T, 

a 



P 


1.00 -0.25 -0.67 


-0.24 


-O.65 


-0.08 


r w -r 


1.00 0.01 


-0.65 


-0.01 

V 


-0.34 


o) 


1.00 


0.45 


0.98 


0.31 


q 




1.00 


0.48 


0.36 


cOq 






1.00 


0.28 
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TABLE V 



Regression equations using only 
the first variable entered. 





Case 


Variable 

entered 


R 


Regression 

coefficients 

A o A 1 


Root mean 
square error 
of estimate 



A 


o>q 


0.51 


0.57 


-1.66 


4.65 


B 


u> q 


0.65 


-0.01 


-0.32 


1.32 


A 


u) 


0.43 


1.07 


-9.20 


4.89 


B 


u) 


0.67 


-0.08 


-1.37 


1.28 



TABLE VI 



Matrix showing the simple correlation 
coefficients between the variables for 
Case A. The primed variables were computed 
for the beginning of the time period. 





P 


Or r 


to 7 


q' 


cb'q' 


T'-T' 


P 


1.00 


0.28 


-0.65 


0.04 


-0.64 


- 0.10 


rv-r 




1.00 


-0.24 


-0.39 


-0. 14 


-0.30 


u>' 






1.00 


-0.40 


0.99 


0.36 


q' 








1.00 


-0.42 


-0.62 


edV 










1.00 


0.33 
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TABLE VII 



Same as Table VI for Case B 





P 


r„-r 


u)' 


q' 




t'-t' 


p 


1.00 


-0.25 


-0.62 


-0.11 


-0.66 


0.37 


i 

I 




1.00 


0.06 


- 0.64 


0.04 


- 0.04 


tX >‘ 






1.00 


0.44 


0.93 


- 0.32 


q' 








1.00 


0 . 4 l 


0.01 


u) / q / 










1.00 


-0.26 



TABLE VIII 

Same as Table V for variables 
computed for the beginning of 
the time period. 


Case 


Variable 

entered 


R 


Regression 

coefficients 

A o A 1 


Root mean 
square error 
of estimate 


A 


U)'q' 


0.64 


1.71 


-2.08 


4.16 


B 


u)' q' 


0.66 


0.21 


-0.72 


1.31 


A 


U )' 


0.65 


2.22 


- 13.77 


4.14 


B 


CO' 


0.62 


0.07 


- 2.59 


1.36 
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TABLE IX 



Improved regression equations using 
only the first variable entered. 



Case 


Variable 

entered 


Regression 

coefficients 

A o A 1 


Root mean 
square error 
of estimate 


A 


u) q 


0.20 


- 1.80 


4.63 


B 


ui q 


O 

i-H 

• 

r-i 

1 


- 0.50 


1.19 


A 


u) 


1.00 


- 0.35 


4.88 


B 


u) 


- 0.80 


- 1.85 


1.20 
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4. CONCLUSIONS AND REMARKS 

Probably the most that could be claimed for the 
results of this study is that it may have some 
applicability over the central and eastern United States 
for the winter months during the period 1200 to 0000 GMT. 
However there is indication that the product of vertical 
velocity and specific humidity at the top of the 
boundary layer can be used to predict convective 
precipitation. 

While this study is diagnostic in nature, it uses 
data that can be predicted numerically. Consequently, 
the technique could be included in a prognostic scheme 
with some expectation of success. The technique of 
using lapse rates, as described here, to determine areas 
of possible convective activity could be applied to 
local forecast areas using forecast temperature fields. 

For a synoptic situation to be used in this type 
of a study, there must be a clear delineation between 
large-scale and convective precipitation. To find out 
if a situation meets this requirement, a prediction of 
precipitation amounts due to large-scale processes must 
be made and compared to the reported precipitation. The 
time required for this procedure limited this study to 
the two cases considered here. Furthermore, the small 
size of the samples virtually precludes a statistical 
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proof of the significance of the results. However, 
the findings of this study are sufficiently encouraging 
to justify further work. 



> 
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